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Abstract: Dirichlet walls - timelike boundaries at finite distance from the bulk on which 
the induced metric is held fixed - have been used to model AdS spacetimes with a finite cut¬ 
off. In the context of gauge/gravity duality, such models are often described as dual to some 
novel UV-cufoff version of a corresponding CFT that maintains local Lorentz invariance. 
We study linearized gravity in the presence of such a wall and find it to differ significantly 
from the seemingly-analogous case of Dirichlet boundary conditions for fields of spins zero 
and one. In particular, using the Kodama-Ishibashi formalism, the boundary condition that 
must be imposed on scalar-sector master field with harmonic time dependence depends ex¬ 
plicitly on their frequency. That this feature first arises for spin-2 appears to be related 
to the second-order nature of the equations of motion. It gives rise to a number of novel 
instabilities, though both global and planar Anti-de Sitter remain (linearly) stable in the 
presence of large-radius Dirichlet cutoffs. The instabilities arise on the outside of spherical 
Dirichlet walls, and also inside sufficiently large such spherical walls in de Sitter space. We 
analyze both inside and outside of flat and spherical walls in Minkowski, de Sitter, and 
anti-de Sitter space, as well as in certain black hole spacetimes and find stability for cases 
not mentioned above. In particular, we find no linear instabilities in the presence of flat 
walls. We also find evidence supporting the conjecture that neutral black holes are repelled 
by Dirichlet walls. 
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1 Introduction 


It is often interesting to consider only a portion of a given spacetime, a construction that 
in field theory is sometimes known as “putting the system in a box.” See [1] for an early 
application of this idea to gravitational physics. The introduction of geometric cutoffs or 
hard walls in a given spacetime is commonly used to model more realistic situations in 
which the wall, or rather, the boundary conditions whose presence the wall requires, is 
the result of some physical property [2-4]. The point of doing so is to typically to make 
a mathematically simpler toy model of a complicated phenomenon. Geometric cutoffs are 
also widely used in the context of gauge/gravity duality, in which they not only play the role 
of dual ultraviolet [5-10] and infra-red cutoffs [11], but also provide toy models of physical 
boundaries in the dual gauge theories [12]. (See [13] for more exact descriptions.) In the 
particular case where the intrinsic metric is fixed, we call this surface a “Dirichlet wall.” 

With such motivations in mind, we recently considered the problem of a black hole 
moving through a spacetime with a Dirichlet wall [14]. In order to make the problem ana¬ 
lytically tractable, we considered the extreme multi-black-hole solutions of [15, 16]. Using 
the method of images to satisfy the boundary condition at the wall, one may obtain ge¬ 
ometries describing static black holes in the presence of a flat Dirichlet wall. Their relative 
motion may then be studied perturbatively in a small velocity (or moduli space) approxi¬ 
mation as in [17]. Though our analysis raised questions as to whether this approximation is 
really well-defined, we found robust evidence that the kinetic energy of such extreme black 
holes can become negative near a Dirichlet wall. This signals a likely dynamical instability. 

However, [14] found no evidence of instability for a purely gravitational theory in 
which all black holes are uncharged. To gain further insight into gravity with Dirich¬ 
let boundary conditions, we now investigate the linear stability of vacuum gravity about 
various Dirichlet-wall backgrounds^. We choose backgrounds constructed from the usual 
maximally-symmetric spacetimes (Minkowski, Anti-de Sitter (AdS), and de Sitter (dS)) in 
arbitrary numbers of dimensions by adding spherical or planar static Dirichlet walls. The 
planar cases turn out to be stable. In the spherical cases we can distinguish perturbations 
outside the wall from those inside, with the latter being confined to a spherical cavity. In 
stark contrast to the familiar results for spins zero and one, the systems outside spherical 
cavities are all unstable. Turing to the situation inside the cavity, the Minkowski and AdS 
systems inside spherical walls are stable though the corresponding dS system is unstable 
for sufficiently large spherical walls. In some situations these linear stability results can be 
established analytically. When this is not the case we resort to numerical computations. 

For the systems inside spherical walls (and inside a flat wall in AdS), we also study 
linear perturbations in the presence of black holes. We find new modes that describe stable 
oscillations of the black hole within spherical cavities in either Minkowski or AdS space, 
thereby confirming the intuition that horizons promote stability. This result also supports 
the prediction of [14] that uncharged black holes are repelled from Dirichlet walls. We also 
find the expected instability to roll down the gravitational potential toward the cosmological 

^To be precise, this manuscript will be entirely devoted to the study of linear-mode stability of Dirichlet- 
wall backgrounds. 
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horizon in de Sitter space - provided the black hole begins far enough from the Dirichlet 
wall that the above repulsion can be neglected at short times. 

Our analysis relies heavily on the Kodama-Ishibashi formalism of [18], which decouples 
the various types of gravitational perturbations (known as scalar, vector, and tensor) by 
writing them in terms of so-called master fields. An important part of our work below is 
to determine how Dirichlet boundary conditions on the metric should be implemented in 
terms of the master fields. The boundary conditions for vector and tensor master fields are 
standard ones (e.g. Dirichlet for the tensor modes), and as a result these modes are stable^. 
But the boundary condition for (quasi-)normal modes of scalar-sector master fields turn 
out to explicitly depend on the frequency. This leads to a very different type of eigenvalue 
problem than is familiar from the study of spin-0 or spin-1 fields and can lead to instabilities 
even when both tensor- and vector-modes are completely stable. As will be described below, 
the fact that such novel boundary conditions first arise for spin-2 appears to be linked to 
the second-order nature of the equation of motion. 

This article is organized as follows. Section 2 reviews the master field formalism of [18] 
and determines the approprate boundary conditions on master fields at a Dirichlet wall. 
We also present an algorithm to reconstruct the full metric perturbation from the master 
field and its derivatives. Finally, we show that adding a black hole to our system leads to 
the existence of new modes that can be obtained by acting with diffeomorphims but which 
are nevertheless physical. These modes are manifestly stable except in those de Sitter 
cases where the gravitational potential is outward-decreasing at the Dirichlet wall. Such 
cases occur when the wall is (sufficiently) closer to the cosmological horizon than to that 
of the black hole. Section 3 presents our results for spherical domain walls in Minkowski 
space, AdS, and dS. In section 4 we consider perturbations in the presence of flat walls in 
pure Minkowski space and in AdS. We summarize our results in section 5, discussing some 
of their implications and the relation to certain non-perturbative issues. A few aspects 
or our treatment, and in particular analyses of non-diffeomorphism modes in black hole 
spacetimes, are relegated to the appendices. 

2 Master field formalism 

The linear stability of gravitational solutions has long been an important subject of inves¬ 
tigation. Early studies date back to [23-25], where four-dimensional gravitational pertur¬ 
bations where decoupled. With the aid of these techniques, the stability of Schwarzschild 
black holes under linear gravitational fluctuations was finally established in [26-28]. Over 
the years, similar formalisms to explicitly analyze linear stability have been developed [29- 
32], all based on reducing the system of perturbations to decoupled scalar equations. 

^The master field formalism has been used to show stability of black hole solutions without Dirichlet 
walls in [19-22]. These works take the fluctuations (and hence the master fields) to have compact support, so 
that they can show linear stability by constructing the appropriate self-adjoint operators whose eigenvalues 
are the squared frequencies of the Fourier modes. The Dirichlet boundary condition satisfied at the wall by 
our tensor-mode master fields often allows us to adapt this argument to show their stability. But for the 
vector and scalar master fields this argument can be obstructed by the fact that do not vanish at the wall. 
We therefore rely on other methods as well. 
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The most general such formalisms for both neutral and electrically charged spacetimes 
are due to Kodama and Ishibashi [18], [19]. These allow one to deal in a gauge invariant way 
with perturbations in arbitrary dimensions, in spacetimes that are warped products over 
any maximally-symmetric space, for all values of the cosmological constant. Making use of 
this decoupling, many black hole spacetimes with different asymptotics have been shown 
to be linearly stable [19-22]. Alternatively, it is sometimes possible to use local Penrose 
inequalities to argue that some instability must exist [33-35] though we will not use this 
method below. 

This section reviews the master field technology of [18] and expresses our Dirichlet 
boundary conditions on metric perturbations in this language. We also reconstruct the 
metric fluctuations in terms of the master fields, making it explicit that for a given solution 
of the master field equation satisfying the appropriate boundary condition there is a lin¬ 
earized solution of Einstein’s equations satisfying our Dirichlet boundary condition at the 
wall. Finally, we study the pure-diffeomorphism modes that arise when places a black hole 
inside a spherical cavity formed by a Dirichlet wall. 


2.1 Review 

Consider an (n -|- 2)-dimensional metric of the form 

= gf.udx^dx’" = gab{y)dy°‘dy^ + r'^{y)dal, 


( 2 . 1 ) 


where gab is a two-dimensional Lorentzian metric and da^ = aijdz^dz^ is the metric on 
the maximally symmetric base space with unit curvature K = 0, -|-1. We leave the case 
K = —1 for future work^. Coordinates x^ labelled by Greek indices are space-time coordi¬ 
nates, which we split them into base coordinates, z*, and orbit coordinates y°‘. We will be 
specifically interested in static solutions of vacuum gravity with metric 


ds^ = —f{r)dt^ + 


dr"^ 2j 2 


f{r) = K - T1 _ Arr 

^ ^ r^n—1 


( 2 . 2 ) 


Here aij is the metric on a n-dimensional unit sphere or an n-dimensional plane. In the 
former case, Einstein’s equations fix X = 1, while in the latter they imply X = 0. The 
constants M and A are related to the total energy E and the cosmological constant A by 


E = 


nMAr. 

SttG 


A = 


n{n + 1)A 


where 


An= dEy^Ja 


(2.3) 


(2.4) 


is the area of the base space. The overall strategy of [18] is to organize the perturbations 
of (2.2) according to how they transform under isometries of a. In practice, this amounts 
to writing a generic perturbation as 


V = 


fJ.U 




(2.5) 


®Black holes with hyperbolic horizons have previously been studied in [36], and their linear stability was 
recently established in [37]. 
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where S, V and T indicate the scalar, vector and tensor parts of the metric, each being a 
sum of terms of the form 


'■ 1 ? = /l?s. 
'> 11 ’= 0 , 
'■ 1 P= 0 . 


'll? = r/f Is,. 

'•S'’ = 

'. 1 T’= 0 , 


hf = 2r>'(J?f Vi,S + fff’Si,), 

/i.S'’ = 2r^H!^h,j, (2.6) 

= 2r24^lTi,-, 


where S, Vj, Tjj denote scalar, vector and tensor harmonics, which are functions on a that 
satisfy 


(ip2 + kl)S = 0, (2.7) 

+ kl)Yi = 0, W = 0 (2.8) 

+ k^)T,j = 0, ViTi = 0, T% = 0. (2.9) 

Here Vi is the covariant derivative associated to a with V^ = PjlD* the corresponding 
Laplacian. The eigenvalues kg, ky and k'^ are non-negative. For K = 0 they form a 
continuous spectrum, while for K = 1 they are discrete and take the form 

kl=esies + n-l), £5 = 0 , 1 , 2 ... ( 2 . 10 ) 

kl =iv{iv+ n-l)-I, = ( 2 . 11 ) 

k^ = ij,(ij, + n-l)-2, It = 1,2,... ( 2 . 12 ) 


As noted in [18], the modes for which = 0 and kg = nK are special and we will 
discuss them separately below. We shall refer to modes with kg{kg — nK) 7 ^ 0 as generic 
modes. For these perturbations, the remaining quantities Sj, Sjj, Vjj appearing in (2.6) can 
be written in terms of the basic harmonics as 

Sj = ——ViS, E>ij = rj^ViVjE) + —cJijS, 'Vij = — V(jyjY (2.13) 

The symmetries imply that perturbations belonging to different sectors decouple from each 
other and can be studied independently. To understand the effect of gauge transformations, 
note that inhnitesimal diffeomorphisms generated by a vector held may be similarly 
decomposed into sums of scalar and vector pieces of the form 

= ci^\t,r)§, = L^^\r,t)r§i (2.14) 

ef)=0, er =L('")(r,t)rW. (2.15) 

There are no tensor diffeomorphisms so is diffeomorphism invariant. To extract 

diffeomorphism-invariant information from the scalar and vector sectors, it is useful to 
introduce the quantities 

= 4®’+ iff? ’ + )(0“)Ga, df’=/?’+20(A, (2.16) 

d''’ = /?’ + )((;ff»ff?’. 
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where 



(2.18) 


and Da is the covariant derivative compatible with the metric gab introduced in (2.1). One 
may then show [18] that (2.16) are left invariant by diffeomorphisms generated by vectors 
of the form (2.14). 

The main result of [18] is that the dynamics of the generic modes of each sector is 
encoded in a single scalar equation of the form 




(2.19) 


for gauge-invariant so-called master fields where Og is the D’Alambertian for the 

metric gab and the index I denotes the scalar, vector or tensor sector. The potential terms 
Vi{r) in the master field equation are different for each sector. The master fields are given 
in terms of the basic gauge invariants defined above by 






rf.n / 2—1 


/)^(^’^/2(J)(V)) _ ^n-lp{V)a^ 


where 


with 


rr/ X + 1) 

H[r) = m-\ - X, 


m = kg — nK, 


X = 


2M 

-n— 1 * 


( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 

(2.24) 


The potential Vj for each sector is 

1 


Vsir) = 


16 ^ 2^2 


{—[n'^(n -|- 2)(n -|- l)‘^x^ — 12n^(n -|- l)(n — 2)mx -t- 4(n — 2)(n — 4)m^]y 


-|- n^(n -|- l)^x^ -|- n(n -|- l)[4(2n^ — 3n -|- 4)m -|- n(n — 2)(n — 4)(n -|- l)K]x^ 

— 12n[(n — 4)m -|- n(n -|- l)(n — 2)K]mx + 16m^ -|- AKn{n + 2)m^}, (2.25) 


^v(r) = -g 


Vrir) = 


, 2 — 2)Ar n(n — 2) , n 

K + kv + -- -Xr^ - 


n{n — 2)K n{n + 2) 



(2.26) 

(2.27) 


where y = Xr^. Due to the lack of diffeomorphism-invariant gravitational fluctuations in 
low dimensions, the tensor master fields vanish identically for n < 2 while the vector and 
scalar master fields vanish identically for n < 1. 

Fourier expanding in the time direction with <h(r, t) = we obtain 


dr{fdr4>’'^^) + f 


(2.28) 


For generic modes, equation (2.28) will be our starting point to study stability. 















We now return to the special scalar modes k‘g{kg — nK) = 0. For /c 5 = 0, S is constant, 
and Si, Sij are not defined. This mode simply shifts the values of M and A, so it need not 
be considered further. For K = 1 the mode k^ = riK is clearly distinct from the mode 
ks = 0, and it turns out to have £5 = 1 . As noted in [18], these modes are generated by 
diffeomorphisms. But they may still contain physical information due to the presence of 
the wall^. We will postpone the analysis of this type of excitation until section 2.2.3. 

2.2 Boundary conditions in terms of Master Fields 

As we now demonstrate, the problem of studying the linearized Einstein equation in the 
presence of a Dirichlet wall is equivalent to solving (2.28) with particular boundary con¬ 
ditions. Examining (2.6), we see that fixing the induced metric on the surface r = is 
equivalent to requiring 

/if = /f^ = 4^^=^r=0, f^^^=HP=0, hP=0, (2.29) 

at r = r£). Note that remains unrestricted by our boundary conditions. 

We now derive the appropriate boundary conditions on the master fields which imple¬ 
ment (2.29) in each sector. We also describe how to recover a metric perturbation that 
vanishes at the Dirichlet wall from the resulting master fields. 


2.2.1 Tensor modes 

For the tensor modes, we see from (2.22) that the master field boundary condition is simply 

= (2.30) 

(T) 

Given a solution of the master field equation, the tensor metric perturbation is just = 
with all other fields in (2.29) set to zero. 


2.2.2 Vector modes 

For vector modes, one may note that the boundary condition (2.29) implies 


F, 


iv), 




rHr 


(V) 


= 0 . 


r=rD 


Equation (2.21) then requires. 


= 0 , 


(2.31) 


(2.32) 


which is the desired boundary condition on vector master fields. 

To reconstruct a metric perturbation we use the single gauge parameter available in 
the vector sector (2.14) to set 

= 0. (2.33) 

^These modes are trivial (pure gauge) in the absence of the wall. Perhaps for this reason, they appear 
not to have been previously studied in detail. 
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We then define 


(2.34) 

(2.35) 


(V) _ kv (y) 

Ilqn -- 

iur 

f(V) ^ piv, _ 

where the are meant to be understood as functions of r, and dr(f)^^\ as specified 

by (2.21). The resulting metric perturbation solves the linearized Einstein equation and 
satisfies our boundary condition. 


2.2.3 Generic scalar modes 

We finally turn to the generic scalar modes. These include all scalar modes with K = 0 
and the K = I scalar modes for which is > 1. Consider the gauge invariant quantities X, 
V, Z given implicitly by 


p(s)t ^ (" - 1 )V - r 


nr 


n—2 


F(S)r^ = 


( 5 ), (n-l)V-X 

nrn-2 

X + Y 


-.n—2 ’ 


J?iS) _ _J 


2nr 


n—2 ■ 


Note that they satisfy the constraint 

+ 2(n - 2)F('^) = 0, 


(2.36) 


(2.37) 


which allows us to invert the relations (2.36) to algebraically express X, Y and Z in terms 
of and F^^\ It follows from (2.20) and (2.36) that 


0(5) = 


inujZ — r{X + Y) 
rn/2-lH 


(2.38) 


Moreover, as shown in [18], using Einstein’s equations we can express X, Y and Z in terms 
of 0(5) and 5^0(5). We refer the reader to [18] for the explicit expressions. 

It will prove convenient to introduce 


■= E(5)*i ± 

f 


A simple calculation then gives 


(2.39) 


F- - ^ f( 5 ) , 2ia;r (5) 1 


9 ( 

UJ 2 


ys) rf ( 5 ) 
Ht j. tij^ . 


(2.40) 


We now see from (2.29) that the Dirichlet boundary condition requires E = 0 at r = r£). 
From (2.36), this means that the scalar master field must satisfy the boundary condition 




= 0 , 


(2.41) 


r=rr, 











Recall that X and Y can be written in terms of r, and dr(j)^^^ ■ Hence, the boundary 
condition (2.41) implies a Robin boundary condition for of the form 

a{r)dr(j)^^'> + P{r)(j)^^'^\r=rD = 0 (2-42) 


where 

m{n — 1)/^ 

a = - - - 

rH 

2 kim-nf-H) m(n-l)(2fc| + n(n-2)/)/ m{n-l)H'p 
^ 2nr‘^H rH"^ 


(2.43) 

(2.44) 


Note that the first term in /3 depends explicitly on the frequency w. This novel feature 
appears to arise from the fact that the dictionary relating scalar perturbations to the (ten¬ 
sor) metric field contains extra derivatives, and from the fact that two radial derivatives 
are related to second time derivatives by the equation of motion. 

R is not yet manifest that the boundary condition (2.41) is truly equivalent to fixing 

the induced metric at r = We now show that this is so by giving an algorithm to 

iS) 

construct a scalar metric perturbation hln/ which solves the linearized Einstein’s equations 
and satisfies our Dirichlet boundary condition. First we must choose a gauge. A convenient 
choice is to set 

= 0. (2.45) 


This exhausts the three functions in the general scalar gauge transformation in (2.14). Now 
let the remaining metric functions be given by 


AS) 

Jtt 

AS) 

J tr 

f(S) 

J rr 


-fE- 



(2.46) 


where in these expressions ^re meant to be taken as functions of X, Y 

and Z according to (2.36). Given a solution of the master field equation (2.28) for the 
scalar potential, we can calculate X, Y, Z using the explicit expressions in [18]. We then 
define our metric functions by (2.45) and (2.46). The resulting metric is a solution to the 
linearized Einstein’s equations which satisfies our Dirichlet boundary condition when the 
master field satisfies (2.42). 


2.3 = 1 modes 

As mentioned above, modes with K = 1, is = 1 must be considered separately since they 

(S) 

yield Sjj = 0 identically, so that the function drops out of the analysis above. As noted 
in [18], all solutions with is = 1 can be generated by diffeomorphims. Here we may allow 
any n > 0. A general is = 1 diffeomorphim of the form (2.14) generates the perturbation 
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with Sij = 0 and, dropping the (S) superscripts for the rest of this section since we consider 
only scalar modes, takes the form 


fr 


frt 


ftt 

ft 


Hl 


rV — L — y/ncr 
r ’ 

Cyf + 2fd^ 

f 

fc't - Ctf - iojfCr 

f 

-{crff + 2iuct), 
iojrL + ^/nct 

1 

r 

\/nL + nfcr 
nr 


(2.47) 

(2.48) 

(2.49) 

(2.50) 

(2.51) 

(2.52) 

(2.53) 


Thus for .^5 = 1, the Dirichlet boundary condition requires 


where 


iuroLirn) + y/nctiro) = 0, 

(2.54) 

^/nL{rD) + nf{rD)cr{rD) = 0, 

(2.55) 

(w^ - a;|g=i)T(r77) = 0, 

(2.56) 

, ,2 firn) 

2rn ■ 

(2.57) 


Consider first the case in which we satisfy (2.56) by setting L{rD) = 0. The mode then 
corresponds to a diffeomorphism that vanishes at r = ro- We expect this perturbation to 
be pure gauge since acts trivially on the boundary (see e.g. [38], or [39] for a well-known 
analogous case). In appendix A, we will see by explicit computation in four dimensions 
that diffeos with L^rm) = 0 are indeed null directions of the symplectic structure and can 
thus be regarded as pure gauge. We emphasize that our analysis both above and below is 
independent of the sign of the wall’s extrinsic curvature and thus applies equally well on 
either the “inside” or the “outside” of our wall. 

On the other hand, one may instead solve (2.56) by imposing u = ±Wfg=i. We can of 
course also satisfy any linear boundary condition at the other end of spacetime by choosing 
the diffeomorphism to vanish in that vicinity. The result is a linearized solution generated 
by a diffeomorphism that acts non-trivially on the wall. Such modes can be physical, and 
for M 7^ 0 one may verify that this is indeed the case by showing that they have non¬ 
zero symplectic inner products. Since it differs qualitatively from the rest of this work, 
we have relegated this calculation to appendix A. There it is performed explicitly for the 
four-dimensional case (n = 2), but the result should be similar in other dimensions. In 
contrast, even modes with ut = ±a;£g=i are pure gauge for M = 0. So for M 7^ 0 they are 
naturally interpreted as describing the motion of the black hole’s center of mass relative to 
the spherical wall. Note that they are stable when the gravitational potential / increase at 
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the wall in the same direction as the radius r, while they are unstable when / decreases in 
this direction. 

For M 7^ 0, A < 0 it is straightforward to show that (2.57) is always positive. So in such 
cases the black hole’s center of mass executes stable oscillations. This confirms expectations 
from [14], which argued that uncharged black holes are repelled from Dirichlet walls. 

But (2.57) can become negative for A > 0 where we find 


UJ 


2 


firp) 

2rp 

1 

2pdphLIs 


(n 


1 ) “ Ph) - “^PDPH 


(2.58) 


for pp = rp/Lds and pp = rp/Lds, and where we require pp < pp < ^ so that the wall is 
timelike (and lies inside the de Sitter horizon). The resulting instability is no surprise for 
Ph PP) where it is clear that the black hole experiences the same instability as geodesics 
in de Sitter which can slide off the central maximum of the effective potential toward the 
de Sitter horizon. And it is reasonable that this instability will become stronger as the edge 
of the black hole approaches the de Sitter horizon (when pp, and thus also pp, approaches 
1). Noting that (2.58) is positive for pp = pp shows that bringing the wall close enough to 
the black hole always removes the instability. This again verifies the repulsive interaction 
between the wall and the black hole horizon. 


3 Spherical walls 

We will first address the case of a spherical wall at r = rp in Minkowski, AdS, and dS. For 
each case, we study separately the fluctuations “inside” and “outside” the Dirichlet wall, i.e. 
we consider the spacetimes for which the metric is (2.2) with M = 0 and K = I, restricted 
to the region r < rp and then to r > rp. In the de Sitter case, we consider only the static 
region inside the de Sitter horizon. Since the modes generated by diffeomorphisms were 
analyzed in section 2.3 and appendix A, we confine ourselves here to modes with is > 

In addition to the boundary condition at rp, we impose regularity at r = 0, outgoing 
boundary conditions at Minkowski r = oo and the de Sitter horizon, and asymptotically 
AdS boundary conditions at AdS r = oo. Thus both the Minkowski and dS cases define 
quasi-normal mode problems rather than normal mode problems. The asymptotically AdS 
boundary condition fixes the the conformal metric at r = oo; for gravitational fluctuations, 
this condition can be obtained by carefully taking the r —)> oo limit of the same boundary 
condition that we impose at our Dirichlet wall. Other boundary conditions appear to lead 
to ghosts [40, 41]. 

3.1 Minkowski space 

We begin with Minkowski space {K = 1, A = 0, M = 0 in (2.2)). For this background, all 
three master field potentials are positive definite and can be written 

Vi = ^[4:ii{ii + n - 1) + n{n - 2)]. (3.1) 
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The wave equation can be solved analytically, but for general modes inside the wall - and 
for tensor modes outside - one may also show stability by extending the argument of [ 20 ] 
and keeping track of the new boundary terms at that arise when integrating by parts. 


3.1.1 Inside the wall 

Consider first the region 0 < r < rjj. In this case we begin by introducing the radial 
coordinate = f f~^dr in (2.28), hnding 

“ (3.2) 


Recall that our boundary conditions are regularity at the origin and, depending on the 
sector, we have (2.30), (2.32), or (2.42) at the Dirichlet wall. Studying the wave equation 
explicitly shows that the two characteristic behaviours at the origin are ~ ^nl 2 +ii 

Thus regularity at r = 0 implies = 0) = 0. Following [20] we 

multiply (3.2) by the complex conjugate of which we denote by and integrate 
both sides. After integrating by parts we obtain 


rf*{ro) 

/r*(0) 




ro rr*(ro) , ^ 

0 -'?'t(0) 


The boundary conditions in all three sectors can be written 

[ai{r, ki)dr(p^^'> + 6 /(r, ki)(py'> - u}‘^c(r, kj)(py'>]= 0 , 


(3.3) 


(3.4) 


where a/(r, kj), bj{r, kj) and ci{r, kj) are functions of the radial coordinate and the radial 
coordinate and the spatial momentum. For the tensor modes, bx = 1, ut = ct = 0, so 
it immediately follows from (3.3) that > 0 for any non-vanishing mode. For the vector 
and scalar modes, we can use the boundary condition at the wall (3.4) in (3.3) to arrive at 


Jr,(o) ai[rD,ki) 

[ \(p^^^‘^dn + (3.5) 

\Jr,{0) ai{rD,ki) J 

For the vectors we have 

av = r, 

so using (3.5) we can conclude that 
turns out that 

= ^5 = ^|^[2A:| + n(n-2)], C 5 = 1, (3.7) 
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by = -, cy = 0, (3.6) 

these modes are all stable. Finally, for the scalars it 


which also leads to > 0. Thus pure Minkowski space is stable against all gravitational 
perturbations that live inside a spherical Dirichlet cavity of finite radius. 







3.1.2 Outside the wall 


Consider now the region r > ro- We impose outgoing boundary conditions at r = oo, 
which selects the profiles that behave as ~ The above argument can again be 

used to show stability of the tensor modes since their boundary term vanished at r = 
and the outgoing boundary condition requires any unstable mode to vanish exponentially 
as r —?• oo. But since r = is now the lower limit of integration, for vector and scalar 
modes the boundary term in (3.3) now yields contributions with the wrong sign. For these 
modes we analyze the wave equation in detail. The radial profiles can then be written 

(3.8) 

where is the Hankel function of the first kind and uj = {n — l)/2 + ij. The factor of 
(7ru;/2)^/^ is inserted for convenience. It is useful to note that the Hankel functions reduce 
to products of exponentials and polynomials for half-integer vj. For example 

iujr 

= i „ „ + 3iujr — 3), n = 2,ij = 2, (3.9) 

i(jjr 

= o o — 15a;r — 15i), n = 4,.^/ = 2, (3.10) 

iLor 

= —i , . (cj'^r^ -I- lOzcj^r^ — 45a;^r^ — l^hiujr -|- 105), n = 6,£/ = 2. (3.11) 

To obtain the desired spectra we must implement the boundary condition at r = nu¬ 
merically. Note that, because depends on uj only through the combination wr, we can 
define Cj = ujr£, and express our solutions in terms of this dimensionless quantity. In this 
form the results are independent of r^). 

It is illustrative to begin with the tensor modes, even though these are already known 
to be stable. The boundary condition is simply 

er(w) ;= (/>('^)(tD) = 0. (3.12) 

Solving this numerically in the complex plane, we find that the tensor spectrum consists of 
a finite number of solutions for each v that increases with z^, all of which are stable, i.e. they 
lie in the lower half plane, see figure 1. Note, however, that some frequencies are complex 
and, from the profiles (3.9)-(3.11), this implies that the solutions diverge exponentially at 
infinity. This is to be expected since, as remarked above, outgoing boundary conditions at 
infinity define a quasi-normal mode problem rather than a normal mode problem. 

Plugging (3.8) in (2.32), we conclude that the vector modes are given by solutions of 

ev{oj) '■= (w) -|- = 0 (3.13) 

Studying (3.13) in the complex plane one finds that the vector modes behave similarly to 
the tensor ones. There are again a finite number of modes for each z^, all of which lie in the 
lower half plane; see figure 1. 
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Re(a)) 


Re(tij) 


(a) 


(b) 


Figure 1. External Minkowski tensor and vector modes are stable: Spectra of frequencies 
for the tensor (left) and vector (right) modes for r > td, K = 1, M = \ = Q, n = A. These are 
obtained by solving eriui) = 0 and ev(a)) = 0 respectively, with ij = 2, 3, 4, 5 (blue circles, yellow 
squares, green diamonds, red triangles). 


We now turn to the more interesting scalar modes. The boundary condition at the wall 
(2.42) reduces to 

es{u}) := (uj) + —?)-(^('^)(cc) = 0. (3.14) 

\ 2n n — 1J 

Here we do find instabilities corresponding to a pair of complex solutions located in the 
upper half-plane, as illustrated in figure 2. 

We observe numerically that for any given (n, is) there is always a single pair of unstable 
modes, and that the number of stable modes grows with i and n. We have checked this for 
n < 12 and 2 < is < 15. The profiles of the unstable modes decay exponentially at infinity, 
so they correspond to normalizable modes which grow exponentially with time. 



- 3 - 2-10 1 2 3 

Re(6j) 

(a) 



(b) 


Figure 2. External Minkowski scalars are unstable: Spectrum of frequencies (left panel) 
for scalar modes with r > rjj, K = 1, M = X = 0, n = 2. These are given by the solutions of 
es(d)) = 0, with is = 2, 3, 4, 5 (blue circles, yellow squares, green diamonds, red triangles). On the 
right panel we plot the imaginary part of the unstable modes as a function of is ■ 
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3.2 AdS 


We next study Dirichlet walls in AdS {K = 1, M = 0, A < 0 in (2.2)). The AdS radius, 
which we will denote by L^dS: is given by LAdS = 


3.2.1 Inside 


The analysis of the interior region 0 < r < r£) proceeds much as in the Minkowski case. 
In particular, the linear stability proof of section 3.1.1 applies directly to both tensor and 
vector modes as the the potentials in (2.27), (2.26) are again manifestly positive for all n 
when A < 0 and their boundary conditions are the same as in Minkowski space. 

The scalars modes must be considered more closely. The potential for the scalar master 
field reduces to 


Vs(r) = ^>^s + nin-2) ^ (n - 2)(n - 4) ^ ^ 

‘^^AdS 

For n 7^ 3 this is positive definite, so in this case we again argue as in 3.1.1. The boundary 
condition (2.42) is of the form (3.4) with 


as = 2n{n - l)r{l + r'^/L\ds? ( 3 - 16 ) 

bs = n(n — 2)(n — 1)(1 + + 2A:|[n(l + r‘^/L\^g) — 1 ] ( 3 . 17 ) 

cs = 2nr‘^ ( 3 . 18 ) 

Because these quantities are all positive, equation ( 3 . 5 ) implies > 0 for n ^ 3. It remains 
to consider the case re = 3 where Vs is not positive definite. For general re, the solution to 
the wave equation that is regular at the origin can be written 




1 

2 


{2 + £s + w)) 2 


1 + ^5 + w); - (re + 1 + 2is)', —p^ 


(3.19) 

with p = T/LAdS^ ^ = ^I^AdS 2^1(0, 6; c; x) a hypergeometric function. The frequencies 
are given by setting re = 3 and imposing the boundary condition (3.4). For general pn = 
rjj/LAdS this needs to be done numerically, but we can obtain an asymptotic expression 
for /9£) —)> 00 by expanding the boundary condition near this point, which yields 


-2 


a{r)dr4>^^'^ + /3{r)4>^^'^\r=rij = 


Vanishing of the leading term requires 


l2T{is + 2) 


r[{2 + is+ u)/2]T[i2 +is-Co)/2] 


+ oip-^) 


(3.20) 


= =t(^5 + 2 + 2p), p = 0,1,2,... (3.21) 

Now, the potential becomes negative only for p^ > 4:is{is-\-2) +3, which is rather large even 
for small values of is- Combined with stability at pu = 00, this suggests that there should be 
no instabilities. We have confirmed this by numerically scanning the range 0.1 < po Ck 20, 
2 < < 15, finding good agreement with (3.21). For example, for re = 2, p^ = 6 we find 

that the first normal modes are ±u) = 4.004,6.032,8.096,10.197. Thus we conclude AdS 
inside a spherical cavity to be linearly stable under all gravitational perturbations. 
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3.2.2 Outside 


We now consider the region r G (r£),oo). It is convenient to introduce the dimensionless 
radial coordinate u = p~^, so the range of interest becomes 0 < u < ud with ud = I/pd- 
The fluctuations must preserve the conformal metric at u = 0, which amounts to imposing 
the Dirichlet-wall boundary conditions (2.30), (2.32), (2.42) in the limit ud —^ 0, see [42, 43] 
for an earlier derivation of this condition in the master held formalism. 

The prohles that satisfy the boundary condition at u = 0 take the form 


</)(^) 


= u^/^-\l + u‘^f/^2Fi 




(tD - 




— £t F 2), 2 + n + 1 + oj); —-—; — 

F +'n + <F)', - ; — , 


(3.22) 

(3.23) 

(3.24) 


The usual plots of where the real and imaginary parts of the r = rjj boundary conditions 
vanish again show that tensor and vector spectra to consist only of real frequencies (with 
an inhnite set for each v). The tensor case constitutes a check on our numerics since, as for 
A = 0, we may prove stability as in section 3.1.1. 

The behavior in the scalar sector is more complicated. For a given mode, stability can 
depend on ri£) in an interesting way; see hgure 3. However, for any given tijj it appears that 
we can hnd quantum numbers which exhibit instability. So linearized gravity is unstable 
outside any spherical Dirichlet wall in global AdS. 
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Figure 3. The imaginary part of the frequency of the unstable scalar mode outside of a spherical 
cavity in global AdS4 for £ = 2,3,4 (blue circles, yellow squares, green diamonds). 


Since the appropriate boundary condition for an asymptotically AdS spacetime is sim¬ 
ply the limit of our Dirichlet boundary condition, it is natural to conjecture that the oscil¬ 
latory behavior seen between a Dirichlet wall and the AdS boundary would be qualitatively 
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similar to that between two Dirichlet walls. This is in fact the case and we find similar 
instabilities in the scalar sector between two spherical walls in Minkowski space with radii 
r_ < r+. The imaginary part of this unstable mode is plotted in figure 4. We see again 
that each mode is periodically unstable, but no unstable region remains when we consider 
all modes. 

0.4 
0.2 

0 

3 0.0 
E 

- 0.2 
-0.4 

2 4 6 8 10 

rjr_ 

Figure 4. The imaginary part of the frequency of the unstable scalar mode between two cavities 
in flat space for .^ = 2,3,4 (blue circles, yellow squares, green diamonds). 



3.3 dS 

Our final background is de Sitter {K = 1, M = 0, A > 0 in (2.2)). The dS radius L^s 
is given by Lds = and we restrict ourselves to cases in which the Dirichlet wall fits 

inside the cosmological horizon, rjy < L^s- This makes the tensor and vector potentials 
positive definite. 


3.3.1 Inside 


We begin with the interior region 0 < r < r£). Solutions that are regular at the origin can 
be written 


^(T) ^ piT+n/2^-^ _ 

^{V) ^ 


^{£t + iCo), + iw + n + 1); ^(2 ^t + n + 1)] 


-(^y + icD + 1), 2 + iCb + n); 2 + n + 1); 


2 (^S + iu) + 2), 2 (^5 + ith + n — 1); 2 ^ + 1)5 


(3.25) 

5 

(3.26) 

(3.27) 


where p = r/L^s and Co = ooL^s- The spectra are again given by imposing the boundary 
conditions (2.30), (2.32), (2.42) at p = po = fjo/Lds- As usual, we can readily adapt the 


17 - 










proof of [20] to forbid both tensor and vector instabilities. On the other hand, we do find 
instabilities in the scalar sector when the wall is close enough to the cosmological horizon. 
These instabilities can be understood by looking at (3.5) more closely. First, note that the 
coefficients characterizing the boundary condition (3.4) can be obtained from (3.16), (3.17), 
(3.18) by the analytic continuation L^s = By doing so, we see that bs is no longer 

positive definite, so our usual stability proof fails. It is natural to expect that the critical 
value rcrit for which bs{rcrit) = 0 to be an upper bound such that there are no instabilities 
for r£) < rcrit- This turns out to be the case. For n = 2, Vcrit is independent of £5, with 
Tcrit/LAdS = l/\/2. For higher values of n, the value of rcrit depends on is in a somewhat 
complicated way, but at large is we have 

TcritlLAdS = + 0{if)- (3.28) 

Figure 5 plots the value of pjj at which the first unstable mode appears as a function 
of is- 
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Figure 5. Unstable internal dS scalar modes: The data points are the value of the location 
of the wall pjj at which the first unstable modes appear as a function of angular momentum is, 
while the dashed line is rcrit/L acIS- Results are shown for n = 2 (left) and n = 3 (right). 


This instability is particularly interesting both because it is the only instability we find 
in the interior of a cavity and because it is connected to the ffuid/gravity analysis of [9, 10]. 
As shown in [9], the limit r/) —)• Lds reduces the dynamics of the gravitational field to that 
of an incompressible fluid on the Dirichlet surface. The stress tensor of this fluid is given 
by the Brown-York expression [39] evaluated on the wall; i.e., by 

tAB = {Kab - Kiab - C'^ab) , (3.29) 

where 7^5 is the induced metric on the Dirichlet wall, Kab is the extrinsic curvature, 
and C is an arbitrary constant that fixes the zero point of the energy. According to this 
correspondence the energy density of the ffuid is given by 

-^0° = - A:7o° - C7o°) = - C), (3.30) 
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where I runs over all of the boundary indices except t. Unlike the explicit examples con¬ 
sidered in [9], in the present case the size of the spheres increase as we move from the 
cutoff surface towards the cosmological horizon, which has the effect of changing the sign of 
the Kj^ term in (3.30). This means that the fluid has a negative energy density or, if C is 
chosen to be sufficiently negative, a term with negative pressure. Given this unconventional 
sign it is plausible that the fluid is unstable, though it would be interesting to develop a 
more complete understanding of the fluid analog of this gravitational instability. 


3.3.2 Outside 


Turning to the exterior modes, we have r G or equivalently, p G {pD, !)• We im¬ 

pose outgoing boundary conditions at the cosmological horizon, so that the prohles behave 
as r\j (1 — p) near p = 1. The corresponding solutions are 


^(T) ^ piT+n/2(^-^ _ 

^{V) ^ p^V+n/2(l _ p2yi6j/2_^p^ 
^{S) ^ pis+n/2^-^ _ p2yiCjl2^p^ 


-{£t - iw), -{n + l + ix - ioj)] I - iuj]l - p^ 

^(1 -I- ^(n -I- t'y - iCb)] I - iCi)]l - 

^{2 + £s - ^{n - 1 + is - - ioj]l - p^ 


(3.31) 

(3.32) 


(3.33) 


Imposing boundary conditions at the wall, we hnd no unstable modes in the tensor and 
vector sectors (as usual, the tensor can be proven to be stable as in section 3.1.2). For the 
scalars, however, we hnd instabilities for all p^ < 1, see hgure 6. 
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Figure 6. Unstable external dS scalar modes: Value of the imaginary part of the unstable 
modes as a function of pu, for = 2, 3,4 (blue circles, yellow squares, green diamonds), with n = 2 
(left) and n = 4 (right). For small pD, the instabilities for each is are a pair of complex modes in 
the upper half plane. As we increase pu, Imw and Rew decrease until these modes merge on the 
positive imaginary axis, after which Imw increases approaching a finite value as po 1. 


4 Results for flat walls 

We now consider static spacetimes with planar Dirichlet walls. Thus we set K = i) m (2.2) 
(with M = 0 as usual). The restriction to static such spacetimes limits our backgrounds to 
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empty Minkowski space and Poincare AdS. Recall that for K = 0 the spectrum of transverse 
momenta k'j is non-negative and continuous. 

In all cases we will find the system to be linearly stable. With the given planar sym¬ 
metry, we find it convenient to simply calculate QNM frequencies analytically rather than 
give an abstract argument for stability. 

4.1 Minkowski space 

We first study fluctuations in half of Minkowski space in the presence of a flat Dirichlet wall. 
Due to the symmetries of the problem, it is more convenient to use a Cartesian coordinate 
system rather than the one in (2.2), so our discussion in this section largely departs from 
the Kodama-Ishibashi formalism. We denote one of the spatial coordinates by 2: and place 
the Dirichlet wall at 2 = 0. To make this choice explicit we write the line element as 

ds'^ = r]ABdy^dy^ -h dz'^, (4.1) 

where includes the time coordinate t and the remaining n spatial coordinates. We 
consider a perturbation upon which we demand hAB\z=o = 0. It is convenient to 
impose the following gauge condition 


V = 0 (4-2) 

where 

■= . (4-3) 

To demonstrate that this gauge choice is consistent with our boundary condition it is 
sufficient to note that an arbitrary perturbation can be put into the gauge (4.2) by a 
linearized diffeomorphism generated by the vector ^ defined as the solution to the initial 
boundary value problem^ 


= 0 • (4-4) 

In this gauge the field equation takes the from = 0. Using translations invariance of 

the background spacetime in the y^ directions we can make the Fourier decomposition 

V = I , (4.5) 

where I runs over all of the y^ except t. 

We will now prove stability by showing that u is real. Inserting (4.5) into the field 
equation = 0 gives 

{dl +u‘^ - kik^)h\^Jiz) = 0 , (4.6) 

®A solution can always be constructed by appropriately extending the problem to an entire Cauchy 
surface of Minkowski space, choosing parity odd initial data, time evolving, and then only keeping the 
solution in the region z > 0. 
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which has the solution 


(z) = sin {k^z) + cos {k^z) k^ := \/a;2 - kjk^ , (4.7) 

where sjl^J, cjl^J are constants and our boundary condition requires that 

Cab - = 0. (4.8) 

We wish to construct asymptotically flat solutions out of a superposition of these modes. 
In order to to do so we must require that the modes be plane wave normalizable. This 
implies that the kj and k^ are real. Therefore 

a;2 > kik^ > 0 . (4.9) 

Furthermore if cj = 0 then the components of h^y are constants and only the are non¬ 
vanishing. This solution is equivalent to a trivial rescaling of the coordinates. Therefore 
Minkowski space with a Dirichlet wall is linearly stable to gravitational fluctuations. 

4.2 AdS 

We now consider the Poincare patch of AdS {K = 0, M = 0 and A < 0 in (2.2)). For 
convenience we introduce the dimensionless radial coordinate p = rjLAdS- We mention that 
scalar fields in the Poincare patch of AdS with various geometrical cutoffs and boundary 
conditions were previously studied in [44]. 

Any mode in Poincare AdS can be classified by its “boundary mass” defined through 
rn^dy I = — kf. We refer to the modes with Tn^^dy i > 0 as time-like, those for which 

rn^dy 7 = 0 as null and those for which rn^dy / < 0 as space-like or tachyonic. In this last 
case is not positive definite, so space-like modes imply instabilities. In order to establish 
stability, we only need to rule out these tachyonic modes and we will thus concentrate on 
these excitations. Unlike Minkowski space where both sides of a flat wall are identical, flat 
walls in Poincare AdS admit a natural notion of inside the wall vs. outside. 


4.2.1 Inside 


We begin with the region 0 < p < pD inside the wall. As mentioned above, it suffices to 
restrict our attention to the tachyonic modes rn^^dy / < 0. These modes decay exponentially 
near the Poincare horizon, so the reasoning leading to (3.5) applies. The boundary condi¬ 
tions at the wall for vectors and tensors are the same as those in Minkowski space, which 
implies > 0, in contradiction with our original assumptions of rn^^dy / < 0. Therefore, 
we conclude that there are no instabilities in the tensor and vector sectors. 

For the scalars, the potential reduces to 


-2 I'ks , (n-2)(n-4) 


^5 = L-/d, + 


''AdS \ .^2 ' ^ 


which is positive definite for n ^ 3. The boundary condition takes the form (3.4) 


as = {n — l)r^ 


n 


— 2)(n — l)r^ -|- ksL\ds^ 


cs = L 


AdS 


(4.10) 

with 

(4.11) 
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These are all positive, so that using (3.5) we arrive at w > 0 which once again rules out the 
space-like modes. Since the potential is not positive for n = 3, we need to resort to a more 
detailed analysis in this case. To simplify the notation, we introduce 5 and 

take Re ks > 0. The tachyonic prohle for the scalar modes is given by 

=p-'/'i^(n-3)/2(W/3)- (4.12) 

where Ku{x) is a Bessel function of the second kind. Replacing this into the boundary 
conditions at the wall (2.42) we hnd 


K(n+i)/2{xs) = ^ (4.13) 

where xg = i^s / Pd- R is well-known (see e.g. [45]) that this equation has no solutions in 
the complex plane with Rexs' > 0. It follows that there are no tachyons in the spectra, 
and that the interior region of Poincare AdS remains linearly stable against gravitational 
perturbations in the presence of a Dirichlet wall. 

4.2.2 Outside 

Consider now p G {pg),oo) and introduce the radial variable u = 1/p. As for K = 1 case, 
the boundary condition at n = 0 hxes the conformal metric, which amounts to imposing 
(2.30), (2.32), (2.42) in the limit ud —t 0. The solutions that satisfy this condition are 


^ ^ ^ '^(n+l)/2 (^6ndj/,T'^) 5 

(4.14) 


(4.15) 

(j)^ = U ^ 3)/2(^6ndj/,5'^)* 

(4.16) 


Once again, the boundary conditions (2.30), (2.32), (2.42) at the wall lead to identical 
results in all three sectors, namely 


J{n+i)/2{xi) = 0, (4.17) 

where xj = rribndyiUD with ud = 1/pd- The zeroes of this equation are real [45] so all 
solutions have > 0; there are no instabilities. 

5 Discussion 

Our work above studied the stability of vacuum gravity solutions under linear perturbations 
in the presence of a Dirichlet wall. We considered flat or spherical walls in Minkowski and 
Anti-de Sitter space, and spherical walls in de Sitter. Stability was analyzed in the regions 
both inside and outside such walls. We also studied the addition of a black hole to the 
region inside the wall. For generic modes the dissipative horizon only increases stability 
(see Appendix B), but as discussed in section 2.3, adding the black hole also leads to new 
“non-generic” modes associated with moving the black hole center of mass relative to the 
wall. The stability of these modes (for A < 0, and for A > 0 when the black hole horizon 
is sufficiently close to the wall) provide evidence that neutral black holes are repelled from 
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Dirichlet walls as predicted by [14], 


Our investigation of generic modes relied heavily on the Kodama-Ishibashi formalism, 
which reduced the problem to decoupled ODEs obeyed by three master fields (tensor, vec¬ 
tor, and scalar). For tensor and vector modes the master field equations define standard 
eigenvalue problems for the (quasi-)normal modes frequencies. In particular, the Dirichlet 
boundary condition on the metric forces the tensor master field to vanish at the wall, while 
it imposes a Robin (or mixed type) condition on the vector master field. As one would 
expect from analogy with spin-0 and spin-1 fields, we found no instability in either of these 
sectors for any of the cases studied. However, the scalar sector satisfied a more general 
boundary condition with explicit dependence on the frequency. This feature gave rise to 
instabilites that would not occur for fields of spins zero or one. 

Our main result was to identify when instabilities actually occur. The examples with 
flat Dirichlet walls were always stable. And inside spherical Dirichlet walls the scalar sector 
was unstable only with a positive cosmological constant. This instability can apparently be 
understood by noting that the coefficient bs in (3.4) becomes negative. Indeed, as £5 —)■ oo 
our plots (figure 5) of the onset of instability for each is appear to approach the locus 
where 65 changes sign. On the other hand, the field outside a spherical cavity was unstable 
in all cases considered. In the AdS case, the stability of individual modes displayed an 
interesting pattern that also appears when one studies the cavity between two concentric 
spherical walls in Minkowski space; see figures 3 and 4. 

It is tempting to associate the unstable modes outside spherical cavities with the fact 
that the Dirichlet wall allows the ADM energy boundary term at infinity to be arbitrarily 
negative. In particular, the boundary conditions defined by any spherical Dirichlet wall 
admit negative-mass (perhaps AdS or dS) Schwarzschild solutions with arbitrary M < 0. 
The wall hides the naked singularity and the solution is smooth in the physical region 
outside the wall. However, the same negative ADM energy would arise in planar AdS 
where the system is stable. Furthermore, we note that the system outside the wall has 
two boundaries - the asymptotic one at infinity and the boundary at the wall itself. As a 
result, time-translations along the wall are not gauge transformations and instead define an 
additional class of global symmetries® associated with a separately-conserved energy H^aii 
given by an appropriate integral of the Brown-York stress tensor (3.29). The structure is 
similar to that of a two-sided wormhole where the ADM Energy is separately conserved in 
each asymptotic region. On (perhaps AdS or dS) Schwarschild, extending the computations 
of [39] to A 7^ 0 and using the sign appropriate to being outside the wall gives 

H^aii = ^rD\/firo) = - Ar^,. (5.1) 

This wall-energy is minimized by taking M sufficiently positive so that rn = while 
Hwall becomes large and positive as M —)■ —00. In particular, conservation of H^all seems to 

^Standard terminology (see e.g. [38]) would refer to them as asymptotic symmetries, but this choice of 
words could cause confusion for a diffeomorphism supported near a Dirichlet wall at finite distance. 
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exclude what would otherwise seem to be a natural scenario in which initially-flat Minkowski 
space outside the wall decays to M < 0 Schwarzschild with positive-energy radiation far 
from the wall. The ultimate fate of our instability thus remains an interesting open question. 

In a recent article [14], we used the small velocity approximation to study the motion 
of a spherical extreme electrically-charged black hole in an asymptotically flat spacetime 
that contains a flat Dirichlet wall. Our results suggested that system to be unstable due to 
the appearance of negative kinetic energies for black holes near the wall. But at least at the 
linear stability level, the present investigation of vacuum solutions found no sign of linear 
instabilities for flat Dirichlet walls. Indeed, as noted above, neutral black holes seem to be 
repelled by Dirichlet walls. It would thus be interesting to use the master field formalism 
of [19] to study the linear stability of charged solutions in the presence of Dirichlet walls. 
It would also be interesting to study these problems at the truly non-linear level. We leave 
such questions for future work. 

The nonlinear stability of Dirichlet-wall backgrounds with respect to spherically sym¬ 
metric spin-0 perturbations has received much attention over the past few years [46, 47]. 
In [46, 47], a nonlinear instability was found for arbitrarily small spherically symmetric 
spin-0 perturbations. The origin of this instability has been attributed to the fact that 
spherically symmetric spin-0 perturbations have a fully resonant spectrum at the linear 
level^ [42, 48-50]. But when these systems are only asymptotically resonant at large w, 
they can be nonlinearly stable [49]. One may check that with zero or negative cosmological 
constant our spin-2 perturbations inside spherical Dirichlet-walls approach resonance even 
more slowly as cu —>■ oo. We therefore expect these Dirichlet-wall systems to be nonlinear ly 
stable as well. 
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A Symplectic structure for Dirichlet boundary conditions 

This appendix studies the symplectic structure for the theory defined with the Dirichlet wall 
in order to determine which diffeomorphism-modes {is = 1) are non-trivial and which are 
pure gauge. We note that the symplectic current receives a contribution from the Gibbons- 
Hawking term which plays a crucial role in its conservation - unless one works in radial 

^This means that there exists an infinite number of tetrads of normal mode frequencies satisfying ±tui ± 
UJ2 it CDs = iLD4. 
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gauge where this contribution vanishes. After constructing an appropriate inner product, 
we compute the norms of the diffeomorphism-modes described in section 2.3. 

We begin by recalling the algorithm of [41] for constructing a conserved symplectic 
structure from a well-dehned variational principle for a held theory in the presence of a 
boundary. See also [51-53] for related treatments of covariant phase spaces which do not 
study such boundaries in detail. 

Denoting the (not necessariy scalar) helds by cj), we assume that the action 

S[cf>] = [ Lo+ [ La (A.l) 

Jm JdM 

has an extremum when some boundary condition b{(f)) = 0 is satished. This b can be any 
local functional of the helds (j). Thus 

6S= [ ndb (A.2) 

JdM 

when the bulk equations of motion are satished, where tti, may be called the momentum 
conjugate to b. As usual, we take dM to be the part of the boundary where boundary 
conditions need to be imposed in order to dehne a phase space. In particular, we neglect 
any terms lying at past or future boundaries of the system. 

Writing the variation of the bulk term as 

6Lq = {eoms)6(p + dOo, (A.3) 

the condition (A.2) requires the pull-back of 9q to dM to satisfy 

9o\dM = T^bdb — 6La + dOa (A.4) 

for some 6a- The total derivative d9a does not contribute to (A.2) since we again neglect 
terms lying at any past or future boundaries. As argued in [41], we obtain a conserved 
symplectic structure by taking the symplectic current to be 

j = jo - dja, (A.5) 

where jo and ja are the symplectic currents associated to the potentials 9o and 9a, i.e. 

jo = d29o[6i(j)] - 6 i9o[S 24'], jd = d29a[di(j)] - di9a[62(j)]. (A.6) 

The key point is that, since the anti-symmetric second variation of La vanishes identically, 
taking the anti-symmeric variation of (A.4) implies that j vanishes when pulled back to 
dM {j\aM = 0) and evaluated on variations satisfying the desired boundary condition (so 
that 5b = 0). Since no symplectic current flows though the boundary, conservation of 
the symplectic structure j-^j for hypersurfaces S having boundaries only on dM follows 
immediately from the fact that the bulk contribution to this current is closed [51-53], 
djo = 0. 

We now apply this construction to gravity with Dirichlet boundary conditions. We 
are interested in space-times for which dM is a time-like surface of constant r, with unit 
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normal n^j,. We assume that they can be globally foliated by constant r surfaces, on which 
we introduce coordinates y*, so that the metric can be written in the form 

ds^ = N^dr"^ + jij{dy^ + N^dr){dy^ + N^dr), (A.7) 

where the induced metric on surfaces of constant r is 'jij and N, are the radial lapse 
and shift functions, respectively. Note that the normal satisfies n^dx^ = Ndr. 

For Dirichlet boundary conditions on dM, a suitable action is given by adding the 
Gibbons-Hawking term to the bulk Einstein-Hilbert term 

S= f ^{R - 2 A) + 2 / (A. 8 ) 

Jm JdM 

where — Ufj^riy is the induced metric at the boundary and K is the trace of the 

extrinsic curvature In this covariant notation, tensors on dM are degen¬ 

erate space-time tensors which vanish when contracted with n^. In particular, note that 
7 = det ■yij and 7 / det = 0 . 


The bulk contribution to the symplectic current is the standard one for Einstein-Hilbert 
gravity, which we take to be given by 

Jeh = - (1 ^ 2 ). (A.9) 

See [54, 55] for other choices of symplectic currents that differ from (A.9) by total derivatives. 
A general on-shell variation of the action (A. 8 ) is of the form 


5S= I + 

JdM 


(A. 10) 


where is the covariant derivative compatible with the conjugate momentum is given 
by 

= -{K^''-Kyf^’^), (A.ll) 

and 

c>^ = -yf^PSgp^n^ (A. 12) 

is tangent to dM so that is well-defined. 

Since our boundary conditions are dy^uloM = 0, the boundary contribution to the 
symplectic potential is 

Oh = -Vi9'%xan^, (A. 13) 

The boundary contribution to the symplectic current is obtained by taking the antisym¬ 
metrized variation 

fd = + (1 ^ 2), (A.14) 

and the total symplectic structure is 


n = 



Jeh - 



(A.15) 
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As noted above, general arguments imply that it is conserved. We have also checked this 
result by direct computation. 

The combination 

{^19,^29) = (A.16) 

defines the desired inner product, where * denotes complex conjugation. We now explicitly 
evaluate this inner product for the diffeomorphism modes discussed in section 2.3. For 
simplicity we consider only the case of four spacetime dimensions (n = 2), though other 
dimensions should behave similarly. 

It is well known that the bulk contribution to the symplectic structure can be ex¬ 
pressed as a total derivative for all diffeomorphims, allowing us to reduce the radial integral 
to boundary terms. We focus on the contribution from the Dirichlet wall at r^i, since we 
can always choose our diffeomorphism to vanish in the neighborhood of any other bound¬ 
ary. Due to the orthogonality of the spherical harmonics, our diffeomorphism-modes can 
have non-trivial inner products only with other is = ^ modes; i.e., only with other pure 
diffeomorphisms. 

Recall from section 2.3 that the frequencies are either purely real or purely imaginary. 
For the real case, which occurs when /'{rs)) > 0, we find 

{hg,529) = Li{rD)Ll2{rD), (A.17) 

where we remind the reader that L is the free function controlling the radial dependence 
of the diffeomorphism. The result is non-vanishing for L{r£)) 7^ 0, so such excitations are 
physical as expected. Conversely, modes with L{r£,) = 0 have vanishing inner product with 
all modes and therefore are pure gauge. 

The imaginary frequencies are similar®. Let h± denote the solutions whose frequencies 
have respectively positive and negative imaginary parts. Reality of the norm requires both 
{h±,h±) = 0 and (/i_|_,/i_) = (/i_,/i+)*. We can verify by explicit calculation that this 
is indeed the case, with (/i_|_,/i_) given by the right hand side of (A.17) (which is now 
imaginary by assumption). For M = 0, the norm vanishes identically even for L{r£)) ^ 0, 
so the matrix of inner products in this subset is null which implies that these modes are 
pure gauge. But for M ^ 0 the matrix of inner products for these two modes can be 
diagonalized. One of the diagonal-basis modes has positive norm while the other norm is 
negative, but both are non-zero for ^ 0. So these modes are again physical. We 

emphasize the key role of the boundary contribution to the symplectic structure (A. 15) in 
arriving at these results. 

B Black Holes and Dirichlet Walls 

As noted in the main text, for modes that exist already in the absence of a black hole, one 
expects the addition of a black hole horizon on the physical side of the Dirichlet wall only 

®A similar example was encountered previously In the context of a massive Klein-Gordon field in AdS 
in [40]. 
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to add dissipation and to make the modes more stable. We now verify this expectation 
for both spherical and (for AdS) planar black holes inside respectively spherical and planar 
Dirichlet walls. The ^5 = 1 modes were already studied in section 2.3, so we consider below 
only £ 5 / 1 . 

We denote the location of the black hole horizon hy r = rn, taken to be the appropriate 
solution of /{vh) = 0. We thus consider modes supported \n < r < vd and impose 
ingoing boundary conditions at the horizon. In terms of the master fields, this condition 
takes the form [56] 

= {r — + 0{r — th)), r = rn (B.l) 

where is a constant and T is the temperature of the horizon given by 

^ = / {rn) ^ 
dvr 

B.l Schwarzschild 

For A = 0 and M > 0 the spacetime (2.2) is Schwarzschild, with its familiar regular event 
horizon at rn = We focus on the region rn < r < between the horizon and 

the wall, imposing ingoing boundary conditions at the horizon: cj) (r - rn) i) 

near r = rjj. As usual, stability of vector and tensor modes may be shown by adapting the 
argument of [20] but we attack the scalar modes more explicitly. 

For non-vanishing M the equation of motion is solved by Heun functions, about which 
rather little is known. We find it most convenient to simply solve numerically for these 
functions and the associated scalar quasi-normal modes (QNM) by discretizing the wave 
equation using a Chebyschev grid and solving the resulting discrete eigenvalue problem. In 
this way, we have determined the first five QNM for 2 < n < 5, 2 < Q < 10, 0.1 < r^/rn < 
10 finding no unstable modes, see figure 7. For each set of quantum numbers the modes 
decay at late times and may thus be said to be more stable than for M = 0 (where the 
modes oscillate but neither grow nor decay). 

B.2 Schwarzschild dS 

Spherical Schwarzschild dS black holes, corresponding toiF = l,M >0,A>0in(2.2), 
can be studied in much the same manner. This geometry exhibits two horizons: the event 
horizon, th, and the cosmological horizon Vcosm, with rn < vd- We restrict ourselves to 
fH < < 'f'cosm and consider modes in the region < r < r^- The temperature of the 

black hole horizon is 

Ingoing boundary conditions at the horizon require ~ (r — 

We solve the for the spectra of QNM numerically as described in section B.l for 2 < n < 5 
in units with = 1. We find no evidence of instabilities in the tensor or vector sectors®, 

^Indeed, stability of the tensor sector can again be shown as in [57]. 
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rolfH ro/rn 

(a) (b) 

Figure 7. Imaginary part of the lowest scalar QNM’s for a Schwarzschild black hole with £ 5 = 2 , 
3, 4, 5 (blue circles, yellow squares, green diamonds, red triangles) for n = 2 (left) and n = 5 
(right). The plots display the following expected behaviors: i) For ^ rjj, modes with larger is 
experience less damping due to the stronger potential barrier that shields them from the horizon, 
ii) For th near ru the potential barrier has comparable effects on the is = 2,3,4, 5 modes. In this 
regime, the growth of positive gradient energy with is makes the is = 5 mode the most stable of 
the ones shown. 


but the generic lowest lying modes in the scalar sector become unstable for sufficiently large 
cavities, see hgure 8 . This is expected from the M = 0 results of section 3.3 which found 
unstable scalar modes but only stable vector and tensor modes. It seems that large black 
holes introduce large dissipation that completely removes the M = 0 scalar instabilities. In 
addition to the unstable generic scalar modes, the dS black holes can suffer from instabil¬ 
ities corresponding to the £3 = 1 modes. These occur in a large region of the parameter 
space {ph,Pd) characterized by /'{po) < 0; see hgure 8 . 


B.3 Spherical Schwarzschild AdS 

Spherical Schwarzschild AdS black holes are described by taking K = 1, M > 0, X < 0 in 
(2.2). The black hole hoirzon lies at the unique positive solution rn of /(rjy) = 0. As for 
A = 0, we consider modes in the region vh < r < ru. Ingoing boundary conditions at the 
horizon read (j) ^ [r — r h )~^^^where 

7sch,AdS = - 1 -b (n -b 1) 

Once again, we solve for the QNM spectrum numerically by discretezing the eigenvalue 
problem given by the master held equation with the corresponding boundary conditions. 
We perform the numerics using the dimensionless radial variable p = r/LAdS-, denote 
Ph = fu/LAdS and pd = 1 ^ 0 /LAdS the locations of the horizon and Dirichlet wall in units 
of the AdS radius. 

The tensor and vector modes can again be shown to be stable by adapting [20]. We 
thus consider only the scalar modes in detail. For 2 < n < 5 we have computed the hrst 
four QNM modes sampling over 0.1 < pn < Pd ^ 10, 2 < < 15 , hnding no instabilities. 

In hgure 9 we plot the frequency of the least stable scalar mode with £3 = 2. We hnd that 


' H 


^AdS 


(B.4) 
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Ph Ph 

(a) (b) 

Figure 8. Stability region of Schwarzschild AdS black holes for n = 2 (left) and n = 5 (right). 
In the shaded region the is = ^ diffeo modes have imaginary frequencies are thus unstable, i.e. 
f {pd) < 0. The red circles correspond to positions in the parameter space in which we have found 
unstable QNM for 2 < < 10. The crosses indicate other points at which we have performed our 

numerical search, and where no instabilities were found. The figures suggest that stability of the 
is = 1 modes also implies stability of all modes with is > 2. 


the imaginary part of the frequency remains negative and approaches the asymptotically 
AdS value in the limit r^/LAds 1- 

B.4 Planar Schwarzschild AdS 

We finally consider fluctuations around a planar Schwarzschild AdS black hole inside a flat 
cavity at r = ro. We restrict to the interior region rn < r < and set Lj^ds = 1- The 
event horizon is then located at ru = We impose ingoing boundary conditions 

at the horizon, which require the solutions to behave as 0 ~ (r — near 

r = th- As usual, adapting the argument of [20] shows the tensor and vector modes to 
be stable, so we need only solve for the QNM spectra of the scalar modes. Our numerical 
method is to discretize the wave equation and solve the resulting linear eigenvalue equation. 
We have obtained numerically the first four scalar QNMs for 2 < n < 5 and a wide range of 
values of the spatial momentum ks = kg/rn and locations of the wall = td/th finding 
that all the modes are stable, see figures 10 and 11 for sample plots. For large values of yu, 
our results coincide with those obtained in the full range of the radial coordinate, computed 
previously in [58] for n = 2. 
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Figure 9. Real (a,c) and imaginary (b,d) parts of the frequency of the least stable scalar mode 
with ^5 = 2 of a Schwarzschild AdS black hole in a spherical cavity. We vary the ratio ro/LAd'S for a 
small (a,b) and a large (c,d) black hole. The black dashed line are the results for AdS Schwarzschild 
without a Dirichlet wall (taken from table 2 of [43]). 



(a) 


(b) 


Figure 10. Real and imaginary part of the first scalar QNM for a planar AdS^ Schwarzschild black 
hole inside a cavity of size j/rj = 3 as a function of k. The lowest lying mode is hydrodynamic. 
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Figure 11. Real and imaginary part of the first k = 2 scalar non-hydrodynamic QNM for a 
planar ^^54 Schwarzschild black hole as a function of yu. The dashed line corresponds to a QNM 
computed for the full spacetime, obtained from table 2 in [58]. 
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